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ABSTRACT 

We present two epochs of near-infrared spectroscopy of the candidate red supergiant 
counterpart to RX J004722.4-252051, a ULX in NGC 253. We measure radial velocities 
of the object and its approximate spectral type by cross-correlating our spectra with 
those of known red supergiants. Our VLT/X-shooter spectrum is best matched by that 
of early M-type super giants, confirming the red supergiant nature of the candidate 
counterpart. The radial velocity of the spectrum, taken on 2014, August 23, is 417 ±4 
km s“^. This is consistent with the radial velocity measured in our spectrum taken 
with Magellan/MMIRS on 2013, June 28, of 410 ±70 km s“^, although the large error 
on the latter implies that a radial velocity shift expected for a black hole of tens of 
Mq can easily be hidden. Using nebular emission lines we find that the radial velocity 
due to the rotation of NGC 253 is 351 ± 4 km s“^ at the position of the ULX. Thus 
the radial velocity of the counterpart confirms that the source is located in NGC 253, 
but also shows an offset with respect to the local bulk motion of the galaxy of 66 ± 
6 km s“^. We argue that the most likely origin for this displacement lies either in 
a SN kick, requiring a system containing a > 50 Mq black hole, and/or in orbital 
radial velocity variations in the ULX binary system, requiring a > 100 Mq black hole. 
We therefore conclude that RX J004722.4-252051 is a strong candidate for a ULX 
containing a massive stellar black hole. 
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1 INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are defined as off- 
nuclear X-ray sources with a luminosity greater than the 
Eddington luminosity of a 10 M© black hole (BH); in prac¬ 
tice a limit for the 0.3 - 10.0 ke V luminosity ^ 10^® erg 
s“^ is used (for a review see e.g. Feng &: Soria||20lT ). The 
two main scenarios to explain these high luminosities are (1) 
super-Eddington accretion on to stellar mass BHs or neutron 
stars ( Begelman||20Q2 ) and (2) sub-Eddington accretion on 
to intermediate mass BHs (IMBHs; Colbert & Mushotzky 

T^. 


* Observations based on ESO programme 093.D-0256 


A growing body of evidence indicates that the group 
of sources we call ULXs may contain objects of both kinds. 
A recent discovery of X-ray pulses from M82 ULX-2 proves 


that the accretor in the system is a neutron star (Bachetti 
et al.|[20T4 ). Yet this ULX sometimes reaches an X-ray lu¬ 
minosity of 1.8 X 10^° erg s“^, which is clearly orders of 


magnitude above its Eddington limit. Liu et al. (2013) use 


optical spectroscopic observations of the donor star of M101 


ULXl and Motch et al. (2014) use optical spectroscopic and 


photometric observations of the donor star of NGC 7793 
P13 to calculate a dynamical limit on the BH masses in 
these systems. In both cases they find the accretor is most 
likely a stellar mass BH. These recent discoveries support 


the suggestion of Gladstone et al. (2009) and others that 
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the peculiar X-ray spectral state, seen in several ULXs with 


high-quality X-ray data (e.g. Bachetti et ah 2013 Sutton 
et al.||2013 Pint ore et al.||2014 ), is an ‘ultr aluminous’ state 


and that many ULXs are stellar mass BHs emitting above 
their Eddington limit. 

However, the brightest ULXs, with luminosities above 
10^^ erg s“^ (also known as hyperluminous X-ray sources 
[HLXs]), are still good candidates to host IMBHs (e g. [Far- 


|rell et ^ 2009 Jonker et al. 2010 Sutton et al. 2012). 


The brightest known ULX and strongest candidate for host¬ 
ing an IMBH is ESO 243-49 HLXl, with an estimated BH 
mass of ~ 20000 M© ( Farrell et al.||2009 ). M82 ULX-1 and 
NGC 2276-3c are two other ULXs that are candidate IMBHs 
( Pasham et al.||2014 Mezcua et al.|2015 ). 

The most reliable, model-independent method to mea¬ 
sure the mass of the accretors in ULXs is through dynamical 
mass measurements. The first attempts to do so (Roberts 
et al.|20rT Liu et al.|2012 ) used optical emission lines from 


the accretion discs in these systems. However, these do not 
show periodic motion. A more promising approach is to mea¬ 
sure the radial velocity curve of the donor star, as has been 


done for several Galactic X-ray binaries (e.g. McGlintock & 
Remillai^|1986 ). 


The optical spectra of many ULX counterparts do not 
show stellar features, as they are dominated by emission 


from the accretion disc (cf. Grise et al. 2012 Tao et al. 
2012 Sutton et al.|[20T4 ). Some ULXs may however have 


red supergiant (RSG) donor stars. Since RSGs are intrinsi¬ 
cally very bright in the near-infrared (NIR) it is possible to 
measure their radial velocities spectroscopically out to ~ 10 
Mpc. An additional advantage is that in the NIR the con¬ 
tribution from the accretion disc is lower than in the optical 
regime. Also, because of the large orbital separations in such 
systems, irradiation of the donor star is not expected to play 
an important role (Gopperwheat et al. |2007| . Therefore, we 
performed a photometric survey of nearby {D < 10 Mpc) 
ULXs in the NIR to search for ULXs with RSG donor stars 
( Heida et al.|2014 ). 

In this paper we present the results of NIR spectro¬ 
scopic follow-up of the ULX with the brightest candidate 
NIR counterpart {Ks = 17.2 ± 0.5) in the sample presented 
Heida et al^ ( [2^ : RX J004722.4-252 051 in NGG 253 


(hereafter J0047), at a distance of 3.4 Mpc ( Dalcanton et al. 
2009|. The c ounterpart is also listed in the catalogue by 


Bailin et al.| ( |2011| ) (source ID 47620), with V = 21.58 


and I = 19.96. Given this U-band magnitude and the Ks- 
band magnitude from Heida et al.| ( |2014| , we derive that 
the V — K = 4.4 iL 0.5 colour is consistent with that of an 
RSG, assuming both the V and K magnitudes do not vary 
strongly in time. 

The ULX has been studied in X-rays by for exam- 
ple|B arnar dlpoTol (NGC 253 ULXl), |Sutton et al.| 

(NGC 253 XMM2 ), |Pintore et al.|p0l^ (NGC 253 X-1) and 


Wik et al. (20141 (source 7). It is not a particularly bright 
ULX, with a peak luminosity in the 0 .3 - 10 keV band of 
2.4 X 10^® erg s~^( Pintore et al.||2014 ), and it has been ob¬ 
served to drop well below the ULX threshold in other ob¬ 
servations. Additionally, it shows unusually high variability 
for a disc-like object ( Sutton et al.|2013 ). The spectrum was 
classified by Sutton et al. (2013) as a ‘broadened disc’, an ul- 


traluminous state associated with sources accreting around 


their Eddington limit (Gladstone et al. 2009). This would 
imply a mass of < 20 M© for the BH in this system. 


2 OBSERVATIONS AND DATA REDUCTION 

We obtained two spectra of J0047 separated by about one 
year: one with the MMT and Magellan Infrared Spectro¬ 
graph (MMIRS; McLeod et al.||2012 ) on the Magellan Glay 
telescope at Las Gampanas Observatory, and the other with 
X-shooter ( [Vernet et al.||2011 ) on the Very Large Telescope 
(VLT) UT3 at Gerro Paranal. 


Magellan/MMIRS 

The MMIRS spectrum (proposal ID 2013A-GFA-3) was 
taken on the night of 2013, June 28, with a 2-pixel (0.4") 
slit, the HK grism and the HK filter. This setup yields a 
spectrum from 1.25 to 2.45 /xm, with an average spectral 
resolution (^) of R 1400. We used an ABBA nodding 
pattern with a dither offset step of 20". We obtained 16 spec¬ 
tra with exposure times of 300 s each, for a total time on 
source of 4800 s. The data were reduced through the GfA 


MMIRS pipeline (v0.9.7RG7-20140301; Ghilingarian et al 
2013). This pipeline subtracts the two exposures of one 
dither pair, performs a flat-held and residual sky correction 
and rectihes the 2D spectra. The spectra are wavelength- 
calibrated using the positions of telluric OH-emission lines; 
the rms of the residuals of the wavelength solution is 0.6 A. 
The dither pairs are then co-added and telluric absorption 
lines are removed using a telluric standard star observed 
close in time to and at similar airmass as the target. These 
corrected 2D spectra still have both a positive and a negative 
trace of the object. We use the starlink program FIGARO 
for the further data reduction steps: use the optimal extrac¬ 


tion algorithm by Horne (1986) to extract the spectrum of 


the positive trace from the pipeline-reduced 2D spectrum, 
the spectrum of the negative trace from the inverted 2D 
spectrum, and then add the two to obtain the final spec¬ 
trum. The nodding pattern unfortunately projected the pos¬ 
itive trace of J0047 on top of the fainter, extended negative 
trace of another region of NGG 253. We use a narrow ex¬ 
traction region for the positive trace to minimize the effects 
of this configuration. We correct the wavelength calibration 
of the MMIRS pipeline (that is done using wavelengths in 
vacuum) to match the calibration of our other spectra (that 
are calibrated using wavelengths in air). 


VLT/X-shooter 

The X-shooter spectrum (programme ID 093.D-0256) was 
taken in service mode on the night of 2014, August 23. X- 
shooter has three spectroscopic arms that operate simulta¬ 
neously, providing spectral coverage from the near UV to 
the near-IR. We used slit widths of 0.8" in the UVB arm, 
0.7" in the VIS arm and 0.6" in the NIR arm, resulting in 
resolutions of 6200, 10600 and 7780, respectively. We used 
X-shooter in nodding mode, with an ABBA nodding pattern 
and a nod throw of 5". The integration times for the UVB, 
VIS and NIR arms were 260, 210 and 285 seconds, respec¬ 
tively. The total exposure times for the three arms are 2600, 
2100 and 2850 seconds. 
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In the NIR arm, the signal-to-noise ratio (S/N) per 
spectral resolution element in the J- and H-hands is ~ 10 in 
regions that are not affected by atmospheric emission lines. 
Because the throughput of X-shooter drops rapidly towards 
the K-band, that part of the spectrum is not useful. 

We process the NIR observations of the ULX counter¬ 
part and a telluric standard star at similar airmass observed 
close in time to the target with the X-shooter pipeline work- 


ffow in the Reflex environment (Freudling et al. 2013). 


The pipeline produces flat-fielded, sky-subtracted, rectified, 
wavelength- and flux-calibrated ID and 2D spectra. The rms 
amplitude of the residuals of the wavelength solution is 0.1 
A. The wavelength calibration is done with arcs that are 
taken during the day, not at the same airmass as the ob¬ 
servations. This can introduce a small offset with respect 
to the real wavelength solution. We check this by comparing 
the positions of sky emission lines to the line list of |Rousselot| 


et al. 


(2000) and find an offset of ~ 0.3 A, that we correct 


for. We then use the SPEXTOOLS-task Xtellcor_GENERAL 
( Vacca et al.||2003 ) to correct the spectrum for telluric ab¬ 
sorption features. 

The iTband spectra of RSGs contain several strong 
absorption edges from CO and lines from neutral metals 
that can be used to derive the effective temperature of the 
star.The J-band also contains absorption lines from neutral 
metals such as Fe, Ti and Mg that can be used to derive 
abundances, but they are not very sensitive to temperature 
( Davies et al.|20f0 ). We have verified that these absorption 
lines are present in the J-band, but in the remainder of this 
paper we focus our analysis on the Jf-band. 

In the data from the UVB and VIS arms the S/N of 
the counterpart is very low (< 1). However, there are spa¬ 
tially extended emission lines visible that appear centred 
on a source at ~ 5" East from the ULX counterpart (not 
detectable in the NIR arm) that is present in one of the 
nod positions (the ‘B’ position). From inspection of archival 
VLT/FORSl R and /-band images taken in the night of 
2004, July 23 we know that there is a blue (U — / = 0; 
[Bailin et al.||20lT ) source that lies partially in the slits at 
that position. The presence of this source makes it impos¬ 
sible to process the data in nodding mode. We follow the 
X-shooter pipeline manual and use the xsh-wkfstare, oca file 
to process the data in stare mode. We extract the emission 
line spectrum at the position of the counterpart from the ex¬ 
posures without the blue source in the slit (the ‘A’ position), 
and at a position close to the blue source from the exposures 
in the ‘B’ nod position. For each exposure, the pipeline pro¬ 
duces bias-subtracted, fiat-fielded, sky-subtracted, rectified, 
wavelength- and flux-calibrated ID and 2D spectra. It also 
produces ID and 2D spectra of the subtracted sky emission. 
For the VIS arm, where there are many strong sky emission 
lines, we use these sky lines to check the wavelength calibra¬ 
tion and find no systematic offset. The rms amplitude of the 
residuals of the wavelength solution for both the VIS and 
UVB spectra is 0.06 A. 


3 ANALYSIS AND RESULTS 
3.1 Cross-correlation analysis 

We are primarily interested in two properties of the candi¬ 
date counterpart: its spectral type, to give us an indication 


of the temperature and (via the observed absolute magni¬ 
tude) the radius of the star, and its radial velocity to confirm 
that it is a donor star in a binary that belongs to NGC 253 
and is not a foreground or background object. Therefore, we 
cross-correlate our X-shooter spectrum with high S/N spec¬ 
tra of known RSGs of different spectral types (ranging from 
K1 to M4.5) and metallicities. 

The set of RSG spectra that we use contains one Galac¬ 
tic RSG, GD-60 3621, whose heliocentric radial velocity is 
— 17.7 ± 0.4 km/s ( [Mermilliod et al.|2008| however, as their 
individual radial velocity measurements of the star show in¬ 
trinsic variations from -16.0 to-19.0 km/s, we adopt a larger 
uncertainty of 3 km/s. All uncertainties quoted in this paper 
are 1-cr errors). Its spectral type is Ml.5 I ( Levesque et al. 
[2005 ). This star was observed with X-shooter as part of the 
X-shooter Spectral Library on the night of 2010, February 3 
(see Ghen et al.|2014 for a description of the library). We re¬ 
trieved the pipeline-reduced spectrum from the ESO phase 3 
archive. The spectrum has been wavelength calibrated, but 
not corrected for telluric absorption. The observations were 
done with the O.O^'^ slit in the NIR-arm, so the resolution of 
the spectrum is the same as for our X-shooter spectrum (ob¬ 
servations for both J0047 and GD-60 3621 were done with 
seeing > 0.6'', so the resolution is set by the slit width). 

The other RSGs in our comparison sample are those 
described by Davies et al. (2013). They observed 19 RSGs 


of a range of spectral types in the LMG and SMG with 
X-shooter. These spectra have been corrected for telluric 
absorption, but their radial velocities are not as accurately 
known as that of GD-60 3621. Since these observations were 
done with a wide (5") slit, they have a lower resolution than 
our X-shooter spectrum (R ^ 5000). 

We use the program MOLLY for the analysis of the data. 
After reading the spectra into MOLLY, we first use hfix to cal¬ 
culate and subsequently apply heliocentric corrections. We 
prepare the spectra for cross-correlation following the rec¬ 
ommendation in the MOLLY user guid^ first normalizing 
with a second-order polynomial fit and then subtracting a 
fourth-order spline fit to the normalized spectra. We resam¬ 
ple the spectra to a common velocity spacing (10.85 km s“^ 
pix“^ for our X-shooter spectrum, 118.4 km s“^ pix“^ for 
our MMIRS spectrum) with vhin and then run xcor to cal¬ 
culate the cross-correlation functions of J0047 with the com¬ 
parison sample, xcor also calculates the radial velocity offset 
from the cross-correlation function by fitting a parabola to 
the three peak pixels. This is a generally robust estimate 
but it underestimates the uncertainty on the radial velocity. 
To get a better estimate of the uncertainty we fit a polyno¬ 
mial plus Gaussian curve to the cross-correlation functions 
using mgfit As a sanity check we also calculate the cross¬ 
correlation functions of GD-60 3621 with the LMG/SMG 
RSGs. The parts of the spectrum of J0047 that are most af¬ 
fected by noise from strong telluric emission lines are masked 
out during the cross-correlation. 

The cross-correlation between the //-band X-shooter 
spectrum of J0047 and GD-60 3621 yields a velocity shift of 
434.4 ± 1.0 km . This uncertainty is the statistical error 
on the determination of the position of the peak of the cross- 


^ http://deneb.astro.Warwick.ac.uk/phsaap/software/ 
molly/html/USER_GUIDE.html 
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Figure 1. Sections of the i^-band X-shooter spectrum of J0047 (upper, blue line) and the X-shooter spectrum of CD-60 3621 (lower, 
red line). Both spectra are normalized and smoothed to the resolution of X-shooter. For presentation purposes only, we interpolated the 
spectrum of J0047 over areas that are strongly affected by telluric emission (grey shaded areas). Error bars are plotted for every fifth 
data point. The spectrum of CD-60 3621 is shifted vertically by —0.5. The J0047 spectrum is shifted by —417 km s“^ and the CD-60 
3621 spectrum is shifted by 17 km s“^. The dashed vertical lines indicate the positions of atomic and molecular lines (left to right, top 
to bottom: Si l (1.59 /xm), CO bandhead (1.62 /rm), CO bandhead (1.66 /xm), Al l triplet (1.672, 1.675, 1.676 /rm), OH (1.69 /rm). Mg l 
(1.71 /rm). 


correlation function. The total uncertainty includes the un¬ 
certainty on the wavelength calibration of X-shooter (2 km 
s“^ for both spectra) and the uncertainty in the radial veloc¬ 
ity of CD-60 3621 (3 km s“^). Quadratically adding all these 
uncertainties and adding the radial velocity of CD-60 3621 


gives a velocity shift of 417zb4 km s ^. The cross-correlations 
of J0047 with the SMC/LMC sample yield consistent results. 

We obtain the highest values of the cross-correlation 
function for CD-60 3621 and LMC-136042 (the latter is 
an M3-type RSG). In addition, following the procedure de- 
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scribed by Origlia et al. (1993), we measure the equivalent 
widths of the CO (1.62 fim) bandhead and the Si I (1.59 
/xm) absorption line using the light command in MOLLY (see 
Figure where parts of the spectra are plotted, shifted to 
zero radial velocity with respect to the heliocentric frame). 
The ratio of these lines is sensitive to the temperature of 
the star. We find a ratio of log(VFi. 62 /bFi. 59 ) = 0.15 ± 0.1, 
which is consistent with an effective temperature in the 
range 3000 - 4000 K [Forster Schreiber|2000| with the caveat 
that our spectrum has a higher resolution than those of 
Forster Schreiber (2000) and Origlia et al.j ( [1993 ), and it 


is not known if this will affect the line ratio]. 


The cross-correlation functions for the SMC RSGs all 
have lower values than those for the LMC RSGs. From this 
we conclude that J0047 is an early M-type RSG, with a 
metallicity closer to that of the LMC and Milky Way RSGs 
than to that of those in the SMC, which have a lower metal¬ 
licity. This is not unexpected, as abundance measurements 
of H II regions in NGC 253 yield values for 12 + log(0/H) 
in the range 8.5 - 9.0 ( Webster Smith]] 1983 ), comparable 
to abundance ratios found in H ii regions in the Milky Way 
(e.g. Balser et al.|2011 ). Values found for H il regions in the 
LMC and SMC are 8.4 ± 0.12 and 7.98 ± 0.09, respectively 
(Pagel et al. 


The Magellan/MMIRS spectrum of J0047, due to its 
much lower resolution and lower S/N, can not be used to 
constrain the spectral type of the RSG. Therefore, we cross¬ 
correlate this spectrum only with CD-60 3621, resampled 
to the velocity resolution of the MMIRS spectrum. We ini¬ 
tially mask only those pixels that are strongly affected by 
residuals from telluric emission lines. This leaves 169 pixels 
(from a total of 463) for the cross-correlation. The peak of 
the cross-correlation function is found at a shift of 3.5 ± 0.2 
pixels, or 410 ± 20 km s“^. As a consistency check we also 
cross-correlate the spectra using only regions where we could 
visually identify absorption lines. This leaves us with 63 pix¬ 
els. The shift found in this way is consistent with the first: 
3.6±0.3 pixels, or 430±30 km s“^. Correcting these shifts for 
the radial velocity of CD-60 3621 gives values of 400±20 and 
410 ± 30 km s“^, respectively. However, the uncertainty on 
the position of the peak of these cross-correlation functions 
is calculated assuming they have a Gaussian profile. This is 
the case for the X-shooter data but not for the MMIRS data. 
Therefore, the error on the radial velocity of the MMIRS 
spectrum is underestimated. 


To gain more insight in the uncertainties on the ra¬ 
dial velocities, we perform a bootstrap analysis. We use the 
MOLLY command hoot to produce 1000 resampled copies of 
the MMIRS and X-shooter spectra. We then cross-correlate 
these with the spectrum of CD-60 3621 using xcor, and fit 
a Gaussian to the resulting distribution of radial velocities. 
The width of this Gaussian is a good measure of the uncer¬ 
tainty on the radial velocity. After correcting for the radial 
velocity of CD-60 3621 we find a value of 416.6 ±1.6 km 
s“^ for the X-shooter spectrum. As we expected, this value 
is consistent with the one we obtained directly from the 
cross-correlation. For the MMIRS spectrum, we find a ra¬ 
dial velocity of 410 ±70 km if we use the 169-pixel mask 
(see Figure]^, and 410±30 km s“^ with the 63-pixel mask. 
We conservatively adopt that the radial velocity of J0047 as 
measured from the MMIRS spectrum is 410 ± 70 km s“^. 


3.2 Nebular emission lines 


Although the UVB and VIS spectra have very weak continua 
(UVB: S/N < 1, VIS: S/N 1) there are several emission 
lines visible in the data. In the UVB data, these are the [O 
II] AA3727,3729, H/3 and [O ill] AA4956, 5007 lines. In the VIS 
spectrum we detect [N ii]A6548, Ho, [N ii] A6583 and the [S 
ii] AA6716,6731 doublet. Inspection of the 2D spectra reveals 
that these lines do not fill the whole slit but are extended 
in the direction of — and peak around — the bright source 
that dominates the unused nod position B (this source is 
located ~ b" East of the RSG, but is not visible in our NIR 
image due to its blue colour; for the 1 and 2D VIS spectra see 
Figure]^. After reading the spectra of the single exposures 
into MOLLY and applying hfix^ we average them using the 
task average. We then measure the radial velocity of the 
emission line region in the UVB and VIS spectra separately 
by using the mgfit function to fit a set of Gaussian curves 
with a common velocity offset to the above-mentioned lines. 
The average radial velocity we find in this way is 351 ±4 km 


From the fits to the spectra we also obtain line ratios 
of log([0 iii]/Hyd), log([N ii]/H(a) and log([S ii]/H(a). These 
ratios are plotted in Figure]^ adapted from |Ho| ( |2008| . Close 
to the blue source the line ratios are fully consistent with a 
H II region, while at the position of the ULX the line ratios 
are more like those seen in Seyferts. 

The He ll A4686 line, which acts as a photon counter 
for photons in the 54 - 200 eV band (Pakull V Angebault 
1986|), has been observed in several ULX nebulae (cf. [Pakull 
fc Mirioni[[2002[ [Kaaret Corbel[[2009[ [Moon et al.[[201l'| 


Gutierrez V Moon 2014). We do not detect it in our X- 
shooter spectrum, with a 2 — a (99.7 % confidence) upper 
limit of 8.8 x 10“^^ erg cm“^ (corresponding to a lumi¬ 
nosity of 1.2 X 10^^ erg s“^). 


3.3 Probability of chance superposition 

J0047 is located close to one of the spiral arms of NGC 253 
(see Figure]^. In such an environment, the probability of 
a chance superposition of an unrelated object may not be 
negligible. Using the IRAF task daofind we search for point 
sources in our VLT/ISAAC As-band image (Heida et al. 
2014) that are equally bright or brighter than the apparent 


magnitude of the candidate RSG {Ks = 17.2 ± 0.5). We 
find that there are 50 such sources in our image, that has a 
total area of 7200 arcsec^. The error circle around the X-ray 
position of the ULX has a rad ius of 1.1^^ or an area of 3.8 
arcsec^ (see Heida et al.[|2014 for the calculation of the size 
of the error circle and details of the ISAAC NIR image). 
This means the probability of finding a source as bright as 
the RSG in the error circle by chance superposition is 2.6%. 
However, this image includes part of a spiral arm of NGC 
253 while the ULX is located just outside this arm, where 
the point source density is lower. Therefore, 2.6% is a conser¬ 
vative upper limit. Due to the inhomogeneous distribution 
of stars, the value for the chance superposition probability 
depends strongly on the part of the image that is selected 
for the calculation. Therefore, it is not possible to find an ac¬ 
curate value for the chance superposition probability. With 
an upper limit of 2.6% we cannot exclude the possibility 
that the RSG and ULX are unrelated. Only a robust mea- 
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Figure 2. The histograms show the distribution of radial velocities obtained by cross-correlating 1000 bootstrapped copies of the 
MMIRS spectrum (red) and the X-shooter spectrum (blue) with CD-60 3621. The systemic velocity of CD-60 3621 of 17 km s“^ is 
already subtracted. The best-fitting Gaussian curves to the data are also plotted (red/blue lines). Inset is a zoom-in of the image around 
415 km s“^ to show the X-shooter data more clearly. 


surement of radial velocity variations will prove whether the 
RSG is orbiting the compact object responsible for the ULX 


emission. 


4 DISCUSSION AND CONCLUSIONS 

We obtained low- and medium-resolution 77-band spectra 
of the counterpart of J0047, a ULX in NGC 253. From a 
cross-correlation with spectra of RSGs in the Milky Way, 
LMG and SMG we conclude that the candidate counterpart 
is most likely an early M-type supergiant. 

In the near-UV and visible parts of our X-shooter spec¬ 
trum we detect several emission lines with line ratios that 
are in between those of H ll regions and Seyferts (see Fig¬ 
ure]^. The 2D profile of the emission lines shows that they 
are extended in the direction, and peak at the position, of 
a bright and blue optical source at ~ b" (corresponding to 
~ 85 pc at the distance of NGG 253) from the RSG. This 
source is not visible in our NIR image but was reported by 


Bailin et al. (2011) to have = —8.1 and U — / = 0, con¬ 


sistent with a late B or early A supergiant. Glose to this blue 
source the ratios of the emission lines are different from the 
ratios found at the position of the ULX, and fully consistent 
with a H II region. We do not detect the He ll A4686 line 
that is indicative of X-ray illuminated nebulae, with a 2 — a 
upper limit of 8.8 x 10“^^ erg cm“^ s“^ (corresponding to 


a luminosity of 1.2 x lO"^^ erg s ). For comparison. Moon 
et al. (2011) report a luminosity in this line of ~ 8.8 xTo^ 


erg s“^for M81 X-6, a ULX with an X-ray luminosity simi¬ 
lar to J0047. We propose that the nebula is excited both by 
young stars and the ULX, tracing a region with recent star 
formation where also the RSG was formed. A detection of - 
or a stronger limit on - the He li A4686 line could (dis)prove 
this scenario. 

We use these lines to measure the projected radial veloc¬ 
ity of NGG 253 at the position of the ULX and find a value of 
351 ±4 km s“^. This is compatible with the valu e of 350zb30 
km s~^ reported by Hlavacek-Larrondo et al. (2011), who 


measured the rotation curve of the galaxy using Ha obser¬ 
vations. The radial velocity of the counterpart, as measured 
from our VLT/X-shooter spectrum taken in August 2014, 
is 417 zb 4 km s“^. This proves that it is not a foreground 
or background object but is located in NGG 253, confirm¬ 
ing its absolute As-band magnitude of— 10.5zb0.5 and its 
identification as an RSG ( Heida et al.||2014 ). This absolute 
magnitude is too bright for asymptotic giant branch (AGB) 
stars, that have absolute A-band magnitudes between —6.4 
and —8.2 (Knapp et al. 2003). This makes J0047 a very 


strong candidate for a ULX with an RSG donor star. The 
ultimate proof for this would be the detection of periodic 
radial velocity variations in the RSG. 

With an effective temperature in the range of 3000 - 
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Figure 3. Part of the VIS spectrum with the [N ii]A6548, Ha, 
[N ii] A6583 complex and the [S ii] AA6716, 6731 doublet. Upper 
figure: the flux-calibrated spectrum extracted at the position of 
the ULX. Lower figure: the 2D image of the spectrum. Here it can 
clearly be seen that the emission lines are extended. The black 
horizontal line indicates the position of the RSG (not visible) on 
the detector. 


3900 K (spectral types MO-3, Tokunaga|2000 ), the expected 
radius of the RSG based on its absolute TTs-band magnitude 
of —10.5 is ~ 600 — 1600 R©. Assuming Roche lobe filling 
RSG with a mass of 10 M©, the orbital period expected for 
such an RSG orbiting a BH is 4.5 - 20 years. The appar¬ 
ent radial velocity amplitude depends on this period as well 
as on the mass ratio and the inclination of the system. For 
instance, for an inclination of 60°, this velocity amplitude 
varies from ~ 20 km s“^ in the case of a stellar mass BH, 
to hundreds km s“^ for an IMBH (see Figure]^. We do not 
detect a significant difference in radial velocity between our 
June 2013 and August 2014 spectra (Av = 10 zb 70 km s“^). 
The large uncertainty is due entirely to the June 2013 low 
resolution, low S/N MMIRS spectrum, and prevents us from 
putting strong limits on the mass of a possible BH compan¬ 
ion. We calculate the probability of detecting a velocity shift 
larger than 80 km s“^ (our 1-cr upper limit on the velocity 
shift) for a range of black hole masses and orbital periods, 
taking into account that the orbital phase is unknown. For 
a system inclination of 60°, we find that if the orbital period 
is 4.5 years, we have a < 50% chance of detecting a velocity 
shift if Mbh < 100 M©. If the orbital period is more than 
10 years we would not detect RV variations even if the RSG 
orbits a 1000 M© BH. 

However, we do measure a radial velocity offset of 66 zb 6 
km s“^ between the RSG and its surroundings (as traced by 
the emission lines). We consider three possible explanations 
for this offset: 

1: The RSG and the X-ray source are unrelated, and 
the RSG is a single runaway star. In this scenario the ULX 
has another, fainter donor star that we do not detect. |E1-| 


find that 7.2% of RSGs have a space velocity > 30 km s“^ 
(their table 4). Of those runaway RSGs, ~ 10% have a space 
velocity > 60 km s“^ (the l-cr lower limit to our velocity off¬ 
set; their figure 2). We combine this with the upper limit of 
2.6% that we find for the chance superposition probability of 
the RSG with the ULX. Using these conservative estimates 
(since we only measure the radial velocity of the RSG and 
the space velocity is likely larger) we find that the proba¬ 
bility of a chance superposition of the ULX with a runaway 
RSG is < 0.02%. 

2 : The RSG is the donor of the ULX and the offset of 
66 zb 6 km s~^ is the systemic velocity of the binary: the 
BH (or neutron star) received a large natal kick (NK) and 
dragged the RSG (or its progenitor) with it (Repetto et al. 


2012| ). Some Galactic BH binaries show high peculiar veloci¬ 
ties (cf. Miller- Jones|2014 ). However, these are all BHs with 
low-mass companions. Instead, Gyg X-1, a binary contain¬ 
ing a stellar mass BH with an ~ 30 M© donor star (which 
is of the same magnitude as the RSG mass), has a peculiar 
velocity of only ~ 20 km/s. The difhculty with this scenario 
is that the probability of kicking a BH/NS plus 10 M© donor 
in a wide orbit to a high systemic velocity, while keeping the 
binary intact, is very small. 

In order to quantify this we simulate a population of 
binaries containing the progenitor of the BH with a 10 M© 
star as companion, calculating the effect of the BH formation 
event on the orbital and kinematical properties of the binary. 
The distribution of initial binary separations is taken to be 
flat between a minimum value amin, corresponding to the 
orbital separation at which either one of the two components 
fills its Roche lobe, and a maximum value of 1000 amin- 
We model the BH formation event assuming a mass equal 
to half the mass of the BH is ejected instantaneously from 
the binary. In the most standard formation scenario, BHs 
are thought to be produced via fallback on to the nascent 
neutron stars. In such a scenario, BHs are then expected to 
receive a NK, as neutron stars do. The NK will be however 
reduced by the ratio Mns/ATbh, simply by conservation of 
linear momentum: Unk, bh = Unk, Ns(ALNs/AfBH), where 
Unk, ns is drawn from a Maxwellian distribution as in |Hobb^ 


et al. (2005). The NK has a random orientation with respect 


dridge et al. (2011) calculated the runaway fraction of RSGs 


in a simulated stellar population; at solar metallicity, they 


to the binary orbital plane. Through the laws of conservation 
of energy and orbital angular momentum, we calculate the 
orbital parameters (a, e) of the binary right after the BH is 
formed. We select among those binaries which stay bound 
in the process, those ones which satisfy two constraints. 

The first constraint is that the systemic velocity ac¬ 
quired by the binary in the BH formation event, i.e. the 
velocity of the new centre of mass with respect to the old 
one, has to be larger than 60 km s“^. The second con¬ 
straint comes from the spectral and photometric properties 
of the companion star. The orbital separation after the BH 
is formed has to allow for mass transfer to happen on the 
nuclear time-scale of the companion when it fills its Roche 
lobe as an RSG. This is done making sure that the Roche 
lobe radius of the companion in the circularized orbit is in 
the interval 600 — 1600 Rq. 

We show the results of this simulation in Figure The 
solid line shows the probability for the binary to stay bound 
in the BH formation event. The dashed line shows the proba¬ 
bility for such a surviving binary to acquire a systemic veloc¬ 
ity larger than 60 km s“^. Both probabilities are shown as a 
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Figure 4. The position and error bars of the line ratios in diagnostic diagrams, indicated by black squares (for the spectrum extracted 
close to the blue source) and green pluses (for the spectrum extracted at the position of the ULX; figure adapted from |Ho|20Q8| >. Close 
to the blue source the line ratios are fully consistent with a H ll region, while at the position of the ULX the line ratios are more like 
those seen in Seyferts, indicating that the ULX might play a role in ionizing the part of the nebula that surrounds it. 



Figure 5. NIR image of the region around J0047. Left: 2MASS image of NGC 253. Right: Our VLT/ISAAC As-band image with the 
X-ray position of the ULX indicated by the black circle (radius 95% confidence). 


function of the BH mass. We also show in the plot the range 
of possible systemic velocities acquired by the surviving bi¬ 
nary; the minimum and maximum value for the systemic 
velocity are calculated at the 99.7% probability. These sim¬ 
ulations were done conservatively assuming an RSG of 10 
Mq; if the donor star is more massive, the probabilities will 
be lower. 


Projection effects lower the probability even further. 
The chance that Vgys is pointed within a degree from our 
line of sight is of the order of 10“^, so the space velocity of 
the RSG with respect to its environment is likely larger than 
what we measure in the radial direction. Bearing in mind 
that BH formation is still very uncertain, we also tested a 
more conservative formation scenario for the BH, in which 
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Figure 6. The radial velocity amplitude for a 10 Mq mass donor 
in orbit around a 10, 100, or 1000 Mq black hole (lower, mid¬ 
dle and upper line, respectively) as a function of orbital period. 
A system inclination of 60° is assumed. Assuming a Roche lobe 
filling RSG, the expected orbital period is between 4.5 and 20 
years. Smaller donor radii give shorter orbital periods; a higher 
donor mass requires a higher BH mass to reach the observed ve¬ 
locity. The black dashed line indicates a velocity of 66 km s“^; 
the hashed area is ruled out based on the necessary radius of the 
RSG. 



Figure 7. Probability fraction for the binary to stay bound in 
the BH formation event (solid line), and probability fraction that 
such a surviving binary acquires a systemic velocity greater then 
60 km s“^ (dashed line), both as a function of the BH mass. 
The shaded area corresponds to the range of possible systemic 
velocities acquired by the binary as a function of the BH mass. 


the BH does not receive any NK at birth, and the peculiar 
velocity at birth is a consequence of the mass ejection only 
(Blaauw kick, Blaauw 1961). In such a case, none of the 
surviving binaries acquire a velocity greater than 60 km/s. 

These simulations show that if the velocity difference 
is due to the systemic velocity only, a BH of > 50 Mq is 
necessary to explain the measured offset. 

3: The RSG is the donor of the ULX and (part of) the 
offset of 66 km s“^ is due to the binary motion of the RSG 
around the BH. It is difficult to explain the offset as the 
systemic velocity of the binary, but the velocity amplitude 
of an RSG orbiting a BH can be hundreds km s“^. This 


requires a BH mass of > 100 Mq for a system inclination 
^ 60° (see Figure |^. 

Such a high BH mass seems at odds with what is ex¬ 
pected based on analysis of the X-ray spectra of this ULX. 
However, Pintore et al. (2014) suggest this ULX may be 
viewed at high inclination, which would mean a lower mass 
BH is sufficient to explain the RV offset (Figure]^ shows the 
RV amplitude assuming an inclination of 60°). If beaming 
plays a role in creating the high X-ray luminosities observed 
in ULXs, a high inclination might explain the relatively low 
X-ray luminosity of this source. The combination of these 
two effects might be enough to reconcile the two mass esti¬ 
mates. 

With the current data we cannot distinguish among 
these scenarios for J0047. New, high quality (X-shooter) 
spectra will allow us to measure radial velocity shifts if they 
are present. This will allow us to decide whether this ULX is 
a good target for dynamical measurements of its BH mass. 
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